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ABSTRACT 



We consider dynamical scales in magnetized GRB outflows, using the solutions to 
the Riemann problem of expanding arbitrarily magnetized outflows (Lyutikov 2010). 
For high ejecta magnetization, the behavior of the forward shock closely resembles 
the so-called thick shell regime of the hydrodynamical expansion. The exception is at 
small radii, where the motion of the forward shock is determined by the dynamics of 
subsonic relativistic outflows. The behaviors of the reverse shock is different in fluid 
and magnetized cases: in the latter case, even for medium magnetization, a ~ 1, the 
reverse shock forms at fairly large distances, and may never form in a wind-type external 
density profile. 



Introduction 



Magnetic fields may play an important dynamical role in the GRB outflows (e.g. Lyutikov 



2006, 2009). They may power the relativistic outflow through Blandford &: Znajek (1977) process 



(e.g. Komissarov 2005 ), and contribute to particle acceleration in the emission regions. In this paper 
we discuss the dynamics of the relativistic, strongly magnetized ejecta. The results are based on an 
exact solution of a one-dimensional Riemann problem of expansion of a cold, strongly magnetized 
into vacuum and into external medium of density p cx (Lyutikov, submitted); they are reviewed in 

In application to GRBs, we assumes that the central engine produces jet with density po and 
magnetization a (a = B^/po; magnetic field is normalized by V4tt), moving with Lorentz factor 
7^ 3> 1. In fact, parameters j w and a are not always independent quantities: at small radii, when 
the motion of the ejecta is subsonic, they should be determined together with the motion of the 
boundary, see {<L3_ In a supersonic regime, relation between 7^ and a depends on the details of 
the flow acceleration (e.g. in conical flows we expect 7^ ~ v^)- For generality, we do not assume 



any relationship between a and 7^. The ejecta is moving into external density p e 
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2. Riemann problem for relativistic expansion of magnetized gas 

2.1. Simple waves and forward shock dynamics 

Let us assume that the jet plasma is moving with velocity f3 w towards the external medium. 
We found (Lyutikov, submitted) exact self-similar solution of relativistic Riemann problem for 
the expansion of cold plasma with density po an d magnetic field Bq (magnetization parameter 
a = Bq/pq; magnetic field is normalized by \/4vr), moving initially with velocity v w towards the 
vacuum interface 

so = 02!W 3 

A 2/3 A l/3 
Or,' 

where the Doppler factors 5 a = + /3 a )/(l — f3 a ) are defined in terms of the plasma veloc- 
ity (3, local Alfven velocity (3a, self-similar parameter rj = z/t, initial wind velocity (3 W and the 
Alfven velocity in the undisturbed plasma /3a,o = + c)- These equations give the velocity (3, 

density p = U\po/o (Ua = (3a/\Ji — P\) and proper magnetic field, B = (p/po)Bo as a function 
of the self-similar variable rj = z/t (expansion of plasma starts at t = 0, z = and proceeds into 
positive direction z > 0). We stress that these solutions are exact, no assumptions about the value 
of the parameter a and velocity v w were made. 

Particularly simple relations are obtained for plasma initially at rest expanding into vacuum 
(3 W = 0, 5/3 = 1 (Lyutikov, submitted). The flow accelerates from rest towards the vacuum interface. 
The bulk of the flow is moving with Lorentz factor 7' ~ cr 1 / 3 . The flow becomes supersonic at rj = 0, 
at which point 7' = (ct/2) 1 / 3 . The vacuum interface moves with Lorentz factor ^' vac = 1 + 2a. In 
the observer frame the vacuum interface is moving with 8^ = 5 A0 5 W , which in the limit a, ^y w S> 1 
this gives 

Ivac = 47^(7 (2) 

As the flow expands, the local magnetization 

decreases. At the sonic point ai oc = (a/2) 2 / 3 . 

If there is an outside medium with density p ex , we may identify two expansion regimes. For 
relativistically strong forward shocks, so that the post-shock pressure is much larger than density, 
the Lorentz factor of the CD is 
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(the last approximation assumes (T>1). For weak forward shocks the velocity of the CD approaches 
the expansion velocity into vacuum j vac , Eq. Q. The transition between the relativistically strong 
and weak shocks occurs for 

o \ 1/3 

3 po x 



For a < (Jcrit, the forward shock is weak. 



2.2. Existence of the reverse shock 



For cold unmagnetized jets the reverse shock always exists; it is weak for j w < -y/po / p cx and 
strong otherwise ( Sari fc Piran|[l995 ) . For magnetized jets the conditions for existence of a reverse 
shock are more complicated (see also Giannios et al. 2008 Mizuno et al.|[2009 ). There are, in fact, 
two somewhat different regimes for the existence of a RS in highly magnetized outflows. First, if 
ejecta is supersonic with respect to the CD (in term of Riemann waves, this transition corresponds 
to the case when the location of the FS coincides with the location of the rarefaction wave), a 
strong RS must forms. Secondly, if the ejecta is subsonic with respect to the CD, but moves with 
velocity higher than the CD, slowing of the ejecta is achieved by a compression wave, which may 
or may not turn into a reverse shock. One dimensional compression waves are always unstable to 



shock formation (Landau Sz Lifshitz 1959). In contrast, multidimensional subsonic outflow need 
not form shocks. So, formally, the condition for reverse shock is 7^ > "fcD, but in the range 
7cd < Jw < 27 M ,- v /cr the RS may not form, if a more complicated flow patters are allowed. In any 
case, the RS shock, even if it exists, is weak in this regime. 

Conditions for strong reverse shock (which implies a highly supersonic flow, with velocity much 



larger than the Alfven velocity in the upstream plasma) were derived by Kennel & Coroniti (1984). 



In the frame of the CD, the reverse shock is moving with ( |Kennel & Coroniti||i~984"l ) 

1 



P 



',2 
RS 



1 



forcr > 1 



(6) 



If lw > 1c D, the reverse shock is weak (if it exists) and one should use a more detailed calculations 
of the dynamics of perpendicular shocks of arbitrary strength. 

Thus, for a > 1, the existence of strong RS requires 7^ > 2jcdV&, which using Eq. ffl gives 



(7) 




while a weak RS may exist for icd < 7w < ^IwyO'- 

f3 




V Pcx 



(8) 
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The Lorentz factor of a strong RS with respect to the contact discontinuity is \/a (assuming 
(j>l). The Lorentz factor of the reverse shock in the frame of stationary external medium is then 



IRS 



11 1CD 



1CD 



/_3_\l/4 
U28/ 



7m / Po_ 

174 ^p ex 

. V4 

PO 



1/4 



(7 



if 7c D » 
if 7c D < v ' T 



(9) 



The two cases in Eq. (|9j) correspond to RS moving in the same direction as the CD, jcd > \/cr, 
and the RS moving in the opposite direction than the CD, jcd < \fo. Condition 7cd = x/& gives 




7u 



In this case the reverse shock is stationary in the frame of the external medium. 



(10) 



Let us summarize the main results. If the the ratio of ejecta density to external density is 
/ = po/pex, then weak RS can form for ^ w > Jcd, which gives j w ~ \fJo\ strong RS forms 
for j w > x/a^/cD, lw > \J fo~ 3 . RS shock propagates in the forward direction for jcd > \/o~, 
lw > \foJ~j . Forwards shock is relativistically weak for jcd > cr, ^ w > a^VvT an< ^ becomes 
non-relativistic for jcd ~ 1) 7w < 



3. Dynamics of magnetized flows in GRBs 



In this section we apply the previous relationships to consider dynamics of magnetized flows in 



GRBs, generalizing discussion of Sari &z, Piranl ( 1995 ) to strongly magnetized flows. We will derive 



main results in a thin shell approximation (not to be confused with a thin shell case, see below), 
assuming that the distances between the forward shock, the contact discontinuity and the reverse 
shock are small. The velocity of the shocks and contact discontinuity are determined from the local 
force balance conditions. More precisely, they are determined by the local solutions to the Riemann 
problem of the decay of the discontinuity of the flow: there is no memory in the flow. Thin shell 
approximation is likely to be applicable, since for reasonable GRB parameters the reverse shock 



never stalls while expansion is relativistic, see discussion after Eq. (20). 



The ejecta flow is taken to expands conically and carrying toroidal magnetic field. We assume 
that the central source operates for time At s = A/c (A is the initial width of the launched shell) 
and produces a wind with magnetization a 3> 1 (magnetization a = B 2 / p is twice the ratio of 
magnetic to particle energy in plasma frame; magnetic field is normalized by v4tt), moving with 
the Lorentz factor j w . For spherical expansion (expansion along conical surfaces) of magnetized 
flows into vacuum, the magnetization parameter a remains constant outside the fast magnetosonic 



surface (Michel 1973 Vlahakis Sz Konigl 2003). As we will see, the above assumption (that the 
central source produces a flow with a given j w and a) is not self-consistent for small radii, where 
the reverse shock does not form. In this case of subsonic expansion, the flow dynamics cannot be 
specified ad hoc: it needs to be determined self-consistently with the motion of the boundaries. 
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The wind luminosity is assumed to be -L; so = E[ so /At s where Ei so is the isotropic equivalent 
energy released by the central source. Luminosity is produced in a form of Poynting and particle 
fluxes 

L = 4vrr 2 7 2 (B 2 + p ) = 47^ 5 2 -^ (11) 

We are interested in the case a > 1. For numerical estimates we will use the typical values for long 
GRBs: Li so = 10 51 erg s _1 , At s = 100 s, E- lso = 10 53 erg, j w = 300. External density is p ex = m p n. 



3.1. Forward shock dynamics 



In case of magnetized ejecta, as well as in the hydrodynamical case (Sari & Piran 1995), the 



important scales in the problem (Sedov scale Is (13), energy scale te (14), reverse shock formation 



scale tn (18), reverse shock crossing scale r& (21) and spreading distance r$ (22)) are related by a 



quantity (Sari & Piran 1995) 




(12) 



In the hydrodynamical case, the parameter £ determines whether the reverse shock and the rar- 
efaction wave reach the whole ejecta before most of the energy is transferred to the forward shock, 
£ > 1, or later, £ < 1. The dynamics of magnetized ejecta generally follows the hydrodynamic thick 
case, though the meaning of some radii change (e.g., in case of strongly magnetized ejecta rjy is 
the scale of RS formation) . 



There is a number of typical radii where dynamics of the outflow changes. There is Sedov 
radius 



s 



1/3 



4 x 10 18 cmn~ 1/3 



(13) 



where the ejecta and the swept-up ISM material become non-relativistic. 



There is radius r^, where the ejecta deposits approximately half of the energy or momentum 
to the external medium. For supersonic flows, which reached terminal Lorentz factor ^y w , equating 



energy in the shocked medium ^p ex c 2 r E to the total energy Ei so , gives (Rees &; Meszaros 1992) 

-Eiso \ Is 



2/3 



9 x 10 16 cmn~ 1/3 



(14) 



te depends exclusively on the total energy of the explosion and not on its form (magnetic or 
baryonic). For radii smaller than te the ejecta's and the forward shocks' Lorentz factors remain 
constant and equal to the initial Lorentz factor 7^. For larger radii the flow enters the self-similar 
Sedov-Blandford-McKee stage, with Lorentz factor decreasing according to 



3/2 



(15) 
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In case of pure baryonic flow, and only in that case, is also the radius when the swept-up 
mass equals the eject a mass divided by ^ w 



M 

"fwPexC^ 



1/3 



1/3 



(16) 



outflow 9^ 



Here Ek = Ei so /(l + a) is the energy associated with bulk motion of matter. Only in the case of 
zero magnetization te equals tm-, since in that case Ei so = Ek = Mq^u,. [j For highly magnetized 

a- 1 ^ « 1. 

The above description of the forward shock dynamics is, in fact, applicable only in the so 
called thin shell case, £ > 1 (see Eq. (12). In this case the reverse shock quickly crosses the ejecta, 
which becomes causally connected so that all of the ejecta interacts with the external medium. 
Alternatively, in the thick shell case, £ < 1 (see Eq. (12), the reverse shock does not have time to 
cross the ejecta before the causally connected shocked part starts to decelerate at smaller radius 
ttv, Eq. (18). The Lorentz factor starts decreasing, but since new material and new momentum is 
being added to the shocked part of the ejecta, the ejecta and the forward shock behave effectively 
as a self-similar shock with energy supply 



7 



L 



1/4 



1 



,3/4 
'■s 



(17) 



Since in the thick shell case the Lorentz factor starts to decelerate earlier than in the thin shell 
case, the rate of energy transfer to the external medium is smaller, so that the ejecta gives most of 
its initial energy to the ISM at larger distances > (Sari &; Piran 1995). 



3.2. Formation and dynamics of the reverse shock 

The weak reverse shock may form at (see Eq. (1ST) 



rN 



3L 



1 



H V 2vrp/5 M c 3 



1 I 



3/2 



10 16 cmn" 1 / 2 



RS becomes strong at rRS,stiong ~ <^fN (see Eq. ([7])). 

If the outside medium is stellar wind, strong RS forms immediately if 

r 2\l/4 / .V \-V4 



3 Lv w i n d(J z 
2^ c 3 M 



1/2 j 1/4 1/4 I 
220<J L 51 <W,8 I 10~8 Mo/yr 



(18) 



(19) 



1 The two radii tm and te were confused by 



Zhang & Kobayashi 



(20051, who "define the deceleration radius using 



Ek alone [] where the fireball collects I/7™ of fireball rest mass". According to Zhang & Kobayashi (2005) "only the 
kinetic energy of the baryonic component (Ek) defines the afterglow level", while magnetic energy is transferred at 
unspecified "later" time. This is incorrect ( Lyutikov||2005 1 . 
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where v w i n d,8 is the velocity of the progenitors wind in thousands kilometers per second. For smaller 
j w , no RS forms ever (for weak shocks, one should put a — > 1). 

The reverse shock could stall (in the observer frame) at (assuming <t>1) 



3L 1 I 



3/2 



rRS,stall 



8p eX C 3 O 



3 x 10 21 cmn" 1/2 o- _1 (20) 



Since rRs,staii is typically larger than Is, RS does not stall during the relativistic expansion phase; 
thus, the thin shell approximation is generally applicable. 

In the unmagnetized case, an important quantity is the radius when the RS crosses the ejecta 



EA 

PexC 2 



V 4 = /| /4 A 1 / 4 = 10 17 cmn- 1 /4 ( 2 i) 



In the magnetized case, is still a good approximation for the RS crossing radius, but with 
two cavities. First, a delayed onset of the RS, see (18), delays the RS crossing moment. Since 
r 7v/ r A = £ 3 ^ 2 , this delay is not important for £ < 1 (the thick shell case, generally applicable to 
the magnetized ejecta). Also, for subsonic outflows (see below), r& is the distance where the back 
of the outflow catches with the CD, see §3.31 



Second, magnetized shell is necessarily expanding, so that the tail part of the flow is moving 
with 7 ~ 7cz>/(2\/o : ). The typical shell spreading distance is 



r s ~ A7 2 = 3 x 10 17 cm (22) 

Since the spreading occurs with Alfven velocity, the tail part of the flow catches with the CD in 
the Blandford-McKee phase at r ta u ~ v^A- We stress that spreading of magnetic shell, unlike of 
the cold baryonic shell, is unavoidable consequence of the high internal pressure (spreading of cold 
baryonic shell requires internal motion). This is the reason why magnetic outflows are similar to 
the thick shell case of the baryonic outflows. 



3.3. Dynamics of subsonic expansion 

Thus, for a given 7^ and a, at distances r < a rarefaction wave is launched into the flow, 
while the flow accelerates to Lorentz factors larger than 7^. This implies that our assumption that 
a flow has a given 7^ and a is not justified at r < at these distances the flow is effectively 
subsonic and its dynamics needs to be solved self-consistently, taking into account interaction with 
the external medium. The subsonic outflow may be considered as a collection of outgoing fast 
magnetosonic waves propagating from the central source, which constantly re-energize the FS. The 
outflow may be separated into two stages, which we will call "early" and "late", depending on 
whether or not most of the fast waves emitted by the central source have caught up with the CD 
and their energy has been given to the circumburst medium. The transition between two stages 
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occurs at the moment, which is similar to the shell crossing radius in the supersonic case, except 
that in the case of subsonic expansion the CD is decelerating all the time, but with different laws 
before and after the transition. 

At the "early" stage the CD is constantly re-energized by the fast-magnetosonic waves prop- 
agating from the central source. The motion of the CD is determined by the luminosity at the 
retarded time t': 

Ln(t') ~ Pcc^cvr 2 (23) 
(this is a condition of pressure balance between the wind and the ram pressure of ISM in the frame 



of the CD). For constant luminosity Eq. (23) gives Eq. (17) for the Lorentz factor of the CD. 



This is exactly the same estimate as for the intermediate scale in £ < 1 supersonic flows, since 



Lq ~ £i so c/A; also this is the same scaling as in the case of relativistic fluid reverse shock (Sari 



1997). We stress that in the limit of strong FS, Eq. ([5]), the FS dynamics is independent of the 



composition of the flow, only the total power is important, Eq. (23). 



The early stage lasts for r < r^, Eq. (21). At larger radii the flow enters the self-similar 



Sedov-Blandford-McKee stage. At this stage, only a fraction of the shell is interacting with the 
external medium, while the newly shocked ejecta material keeps adding energy and momentum to 
the shocked shell and the ISM, which evolve, effectively, as a flow with energy supply. Finally, for 



r > ta the shock enters Blandford-McKee stage, with Lorentz factor given by Eq. (15). 



4. Discussion 



In this paper we discuss the dynamics of strongly magnetized outflows in GRBs. We find 
that the evolution of the forward shock driven by strongly magnetized outflows are qualitatively 
the same as in the case of fluid shocks. The definitions of radii rjv, rE and involve only the 
total energy of the ejecta, it's thickness and initial Lorentz factor, and not the information about 



it's content, e.g., parameter a. The typical radii (12) are the same for two flows (c/. Eq. (12) of 



the present paper and Eq. (9-10) of Sari <fe Piran 1995). These similarities may be understood, 
first, by noting that jump conditions in perpendicular magnetized shocks may be reduced to fluid 
shock jump conditions, with an appropriate choice of the equation of state, and, second, by the fact 
that the thin shell approximation is applicable in our case (so that the global conservation of the 



toroidal magnetic flux, which modifies the global flow dynamics (Kennel Sz Coroniti 1984), is not 



important). Another reason for this similarity is that magnetic field behaves in many respects as a 
fluid with internal pressure. The only difference in the dynamics of the forward shocks driven by 
magnetized and fluid flows occurs for supersonic flows, 7^ > \/ajcD, at very early stages r < rjy 
or r < te-i see Fig. Q. Qualitatively, magnetized outflows are similar to thick shell hydrodynamic 
outflow, £ < 1 at r > r^. 

Only at very early times, at r < rjv, the forward shock bears information about anergy 
content: forward shock is coasting with j w =const in the fluid case and decelerating 7 oc r 



-1/2 



m 
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Fig. 1. — Evolution of the Lorentz factor of the forward shock for matter-dominated (Left panel) and 
Poynting flux-dominated models (Right panel). Matter-dominated ejecta coasts with the injection 
Lorentz factor *y w until either te (for £ > 1, thin shell case) or until rj\r (for £ < 1, thick shell case). 
At r^r reverse shock becomes strong. At large radii (r > te or r > r^) the outflow enters the 
Sedov-Blandford-McKee regime. For highly magnetized ejecta the Lorentz factor of the CD and 
the FS initially decreases Tjsm oc r~ 1//2 , changing to Sedov-Blandford-McKee regime Tjsm °c r~ 3 / 2 
approximately at r/s.- Reverse shock is launched at and becomes strong at r^a. Due to internal 
expansion of the magnetized shell, the back of the shell catches with the CD at distance ~ r^\fo. 
This is the reason why at distances close to ta the Lorentz factor starts decreasing below the r -1 / 2 
law. 



the magnetized case. Dynamics of the reverse shock is quite different in case of high magnetization. 



First, the reverse shock forms at a finite distance from the source (Eq. 18), and may not form at 



all in a wind environment, (Eq. 19). This fact may be related to observed paucity of optical flashes 



in the Swift era (Gomboc et al. 2009). (The standard model had a clear prediction, of a bright 



optical flare with a definite decay properties (Sari et al. 1996 Meszaros & Rees 1997). Though a 



flare closely resembling the predictions was indeed observed (GRB990123, |Meszaros Rees||1999l >, 
this was an exception.) 

In addition, at distances < r < crr^, where the RS is weak, the formation of the RS shock 
depends on the details of the flow: RS forms if the flow is strictly radial, but need not to form if the 
the flow pattern is more complicated. We suggest that optical variability often seen in GRBs {e.g. 
GRB021004 and most notoriously GRB080916C) is a reflection of the non-trivial flow patterns and 
the corresponding non-steady RS formation. Also, a recent detection of high polarization in optical 



(Steele et al. 2009) indicates a presence of an ordered magnetic field in the ejecta. 



I am greatly thankful to Dimitros Gianios, Sergey Komisarov and Alexandre Tchekhovskoy. 
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